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We propose a nonlocal scalar-tensor model of gravity with pseudodifferential operators inspired 
by the effective action of p-adic string and string field theory on flat spacetime. An infinite number 
of derivatives act both on the metric and scalar field sector. The system is localized via the diffusion 
equation approach and its cosmology is studied. We find several exact dynamical solutions, also 
in the presence of a barotropic fluid, which are stationary in the diffusion flow. In particular, and 
contrary to standard general relativity, there exist solutions with exponential and power-law scale 
factor also in an open universe, as well as solutions with sudden future singularities or a bounce. Also, 
from the point of view of quantum field theory, spontaneous symmetry breaking can be naturally 
realized in the class of actions we consider. 



PACS numbers: 98.80.Cq, 11.25.Sq, 11.25.Wx 

I. INTRODUCTION 

Many proposals for modified gravity have been invoked 
in the hope of finding new insights into the open issues 
of the standard cosmological model. Among them, theo- 
ries with pseudodifferential operators have been favored 
with particular attention. A reason is that nonlocal the- 
ories can have very different ultraviolet properties with 
respect to ordinary second- or higher-order actions (in- 
cluding popular Gauss-Bonnet extensions) and, hence, 
could play a role near the big bang and as spacetime ef- 
fective formulations of nonperturbative quantum gravity. 
String field theory (SFT) is a concrete realization of this 
notion where pseudodifferential operators of the form 

(1) 

decorate the effective target action of the fields, where 
r* is a constant and □ is the spacetime d'Alembertian. 
The imprint of nonlocal dynamics in the history of the 
early universe, or even as dark energy models, has mo- 
tivated the study of cosmological models inspired by 
open SFT [1-19], the p-adic string [14, 15, 20-24], or 
other nonlocal effective actions featuring the operators 
(1) [25, 26] or inverse powers of the d'Alembertian [27- 
35]. When nonlocality is of the type (1), it can be con- 
veniently manipulated with the diffusion equation ap- 
proach, which has been developed and employed, in ana- 
lytic and numerical fashion, under different formulations 
[10-12, 14, 15, 18, 23, 36-43]. 

A rather common assumption in the literature of non- 
local fields in cosmology is that nonlocality is confined 
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only to one sector of the model, while the others are lo- 
cal. In the case of SFT-motivated actions, the nonlocal 
sector is matter (a scalar field) and gravity is local and 
with Einstein-Hilbert action: 

S = Sr^onioci^) + d^'xV^ R , (2) 

where D is the topological dimension of spacetime, g is 
the determinant of the metric g^^, ^ = 0,. — 1, 
= 8ttG is Newton's constant, and R is the Ricci curva- 
ture scalar. This ansatz has been dictated mainly by the 
urgency of understanding, in the broadest sense, (i) the 
dynamics of the yet-unclear nonlocal scalar field theories, 

(ii) the combined effect of curvature and nonlocality, and 

(iii) its possible consequences for phenomenology, in par- 
ticular, in relation to cosmology (infiation, dark energy) 
and the modification of flat open SFT solutions (can cos- 
mological friction damp the wild oscillations of the OSFT 
solution with marginal deformations? [4, 14, 16, 18, 37]). 

Now that robust analytical and numerical methods 
have been established to solve nonlocal equations of mo- 
tion, it would be highly desirable to address the concep- 
tual inconsistency subjacent to Eq. (2). Not only would 
we like to define a model with nonlocality implemented 
in all sectors (and reproducing standard general relativ- 
ity in the limit of weak nonlocality), but we want also to 
find nontrivial cosmological solutions. Such is the twofold 
objective of this paper. 

The problem of nonlocal gravity can be faced under 
three independent perspectives, one motivated by string 
field theory, one purely phenomenological, and another a 
hybrid approach. In the first case, the Einstein-Hilbert 
action in Eq. (2) is introduced by hand as an educated 
guess on "how the effective SFT action of tachyon might 
look like in the presence of gravity." Gravity is minimally 
coupled with a tachyon-type or p-adic scalar field whose 
action is dictated or inspired by concrete Minkowski cal- 
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dilations. Obviously, a fully consistent effective tachy- 
onic action should be derived from first principles in all 
its sectors. As far as gravity is concerned, the natural 
framework is closed SFT [44-58] , which features the same 
nonlocal operator (1) of open SFT. The subject is rather 
intricate and, unfortunately, effective gravitational non- 
local actions are known only at the linear level [47, 58].^ 

Instead of facing the rigors of closed SFT some toy 
models have been considered, in particular the open- 
closed p-adic tachyonic action [61-63] and a closed SFT- 
inspired tachyon-tachyon model [64]. However, the gravi- 
ton is not included, thus leaving the (third) possibility to 
consider phenomenological actions where the matter sec- 
tor is as close as possible to SFT or the p-adic string [65] . 
while a nonlocal gravitational sector is built from rea- 
sonable requirements (mainly, that it contains the same 
type of pseudodifferential operators as the matter sector). 
Here we shall follow the third path. 

Nonlocal gravity sectors have been constructed with 
inverse powers of the □ operator [27-35] or more gen- 
eral kinetic functions [26], while keeping matter local. In 
[25] a nonlocal total action has been proposed with non- 
minimal coupling between gravity and a scalar field, but 
the dynamical analysis therein docs not go beyond cos- 
mological solutions when the matter sector is switched 
off. 

In Sec. II we adopt the diffusion equation method to 
infer the form of a solvable scalar-tensor nonlocal action 
with pseudodifferential operators of exponential type, 
Eq. (1). This approach is chosen by virtue of its non- 
perturbative character, which does not require truncating 
the theory in order to find solutions, exact or asymptotic. 
Exact nonvacuum solutions of the equations of motion of 
a p-adic-like system will be found in Sec. Ill for cosmo- 
logical backgrounds; their classical stability is checked. 
The exact solutions are stationary along the diffusion flow 
(i.e., the diffusion equation is trivially satisfied) but the 
scalar and Hubble profiles as well as their dynamics are 
nontrivial. Notably, there exist solutions with exponen- 
tial and power-law scale factor when intrinsic curvature 
is negative (Sec. IIIB), as well as a most general class 
of explicit solutions for actions with conformal operators 
(Sec. IIIC). We make some general remarks on nonsta- 
tionary asymptotic solutions in Sec. HID, showing that, 
for natural choices of the potential, the system realizes 
spontaneous symmetry breaking. 

We are interested in a scalar-tensor theory where non- 
local effects are dominant, and such a system finds a 
natural application only in the early, inflationary, or 
bouncing universe (here we ignore applications of non- 
local models to dark energy). For this reason, in most of 
the paper we do not assume the presence of any cosmo- 



^ On the other hand, the local low-energy effective field theory of 
the closed string tachyon-dilaton-graviton system is well under- 
stood also in its cosmological properties [59, 60]. 



logical fluid such as dust or radiation, but nevertheless 
inclusion of extra matter components is briefly discussed 
in Sec. HIE. In Sec. IV the analysis is extended to an- 
other action with a kinetic operator similar to the one 
of the SFT tachyon; de Sitter and power-law exact solu- 
tions are found. Section V is devoted to discussion. In 
the Appendix we recall flat and curved exponential and 
power-law solutions in standard general relativity. 



II. DEFINITION AND DYNAMICS 

A. p-adic-like action 

Let us begin for simplicity with the p-adic action in 
Minkowski spacetime. 



(3) 



where V is the fleld potential, (p has dimension [(p] = D/2 
in momentum units, and r* is a constant of dimension 
[r*] = —2. Quite often one works in units where the char- 
acteristic energy scale M of the system is absorbed (cor- 
responding to a' = 1 units in string theory) and treats 
r* ~ as dimensionless. In the diffusion equation 

method [40, 43], one introduces an auxiliary direction r 
and a scalar field 4>{r, x) which, for the moment, is not 
related to the solutions (f){x) of the original nonlocal ac- 
tion. By definition, this scalar diffuses according to the 
equation 



{n + dr)(j>{r,x)=Q, 



(4) 



where x = and □ = d^d^^ is the usual d'Alembertian 
operator (later generalized to the curved case). Since 

[n,dr] = o, (5) 

the action of the operators (1) is simply a translation 
along r: 



e''* ^(t){r,x) (j){r, x) = cpir - n,x) . (6) 

At this point the system (3) can be solved by a well- 
established algorithm which, however, we will not need 
because we shall concentrate on "stationary" exact so- 
lutions. Anyway, we briefly summarize it here for the 
interested reader. The idea is to choose an appropri- 
ate initial (in r) field configuration (/){0, x) such that the 
scalar evaluated at the end of the diffusion flow at some 
r = n. 



^{ri,x) = ^{x) , 



(7) 



is indeed the scalar fleld in the action and a solution to 
the equations of motion. To check this, one looks at the 
spacetime asymptotic behaviour of the equations of mo- 
tion and verifles if there exist certain values of the free 
parameter ri such that the proflle (7) is a solution. The 
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existence of this solution is not always guaranteed be- 
cause the diffusion equation does not encode information 
on the potential and, in most applications, the result is 
approximate. In practice, however, the level of the ap- 
proximation [depending both on the chosen configuration 
(^(0, x) and potential V] is more than acceptable and the 
found profile (7) is a genuine nonperturbative solution. 
In string theory this is perhaps not surprising [43] . As a 
consequence of Eq. (6), the system becomes localized in 
spacetime and its infinite number of degrees of freedom 
are transferred into the specification of the field config- 
uration 0(0, x). The resulting model has a well-defined 
restatement of the Cauchy problem [40] (this point, in 
diffusing and nondiffusing models with exponential oper- 
ators as well as in other nonlocal models, has been dis- 
cussed also in [15, 24, 41, 66]). Concrete examples are 
the exact solution for the p-adic string [23] and the ap- 
proximate solutions for open string and superstring field 
theory [39, 42, 43]. 

When introducing a nontrivial metric g^j^, we have to 
include also a nonlocal action for it. As far as the diffu- 
sion method is concerned, so far this has not been done, 
and discussions have been limited to systems of the form 
(2) [10] or with more general local gravitational actions 
[18]. The reason has been mainly technical. Suppose 
decorating the gravitational action with the exponential 
operators (1). If we hope to solve the system analyti- 
cally or semianalytically, it is reasonable to expect that 
a metric-derived field will obey the diffusion equation. 
For instance, this could be the metric itself g^j^y{r,x) or 
one of the Ricmann invariants R{g) = i?^j/...(r, x); dif- 
ferent choices are physically inequivalent. Whatever the 
choice, the fundamental commutation property (5) is no 
longer valid, because now the d'Alembertian is coordi- 
nate dependent and the diffusion equation for the metric 
or Riemann invariant R{g) is essentially nonlinear: 

VA9W^g)R{g) = -drRig), (8) 

where Vo- is the covariant derivative. Accordingly, 

D^R{g) = -DdrRig) ^ -drURig) 

n^R{g) ^ i-irdl^Rig), 

and the translation property (6) breaks down for R{g). 

However, our main goal is not really to impose the dif- 
fusion equation in all sectors, but rather (i) to include 
an infinite number of derivatives of the metric field for 
consistency with the matter sector, and (ii) do it in while 
preserving the translation property wherever and when- 
ever required. Therefore, we can look towards another 
direction: namely, giving up diffusion in metric fields. 



drR{g) = , 



(9) 



and modifying the diffusion equation with an extra curva- 
ture term. Here we take for simplicity a term depending 



only on the Ricci scalar: 

[n + f{R) + dr](l> = 0., 



(10) 



where / is a function of dimension [/] = 2. Equation (3) 
is modified as 



/ 



[□+/(«)] 



(11) 



The role of the ,f(R) term is to realize nonlocality in 
the gravity sector: an infinite number of derivatives act 
on metric fields via the Baker-Campbell-Hausdorff for- 
mula. This way the same exponential operator acts on 
both sectors and one can still solve the equations of mo- 
tion nonperturbatively and analytically, either exactly or 
approximately. 

Equation (11) is the total action for our scalar-tensor 
model. The scalar sector difiiises while the metric sector 
does not, but both sectors are nonlocal as desired. The 
resulting equations of motion can be localized [(j) trans- 
lates, e'''['-'+/(^)l(/)(r, a;) = (f){r — r*,a;)] and solved for a 
given fixed metric. 

Equation (11) has well-defined limits as a pure scalar 
field or gravitational theory. Let /(O) = 0. In the 
Minkowski limit, Eq. (11) reduces to Eq. (3). On the 
other hand, when <p relaxes to a constant (po [local mini- 
mum of the potential, Vq = V{(j)o) = (/)qA] one has, up to 
an overall constant, 



,r,[a+f{R)] 



2A 



(12) 



which becomes an ordinary f{R) theory in the small r* 
limit. In the same limit but keeping (j) dynamical, 

5~ yd^a;y^|i<^[r.n + r./(i?)](/.-y(<^)| . (13) 

Later on we will find solutions of the full nonlocal action 

for the linear case f{R) = aR, where a is a dimension- 
less constant which will be often set to a negative value. 
To recover the Einstein- Hilbert Lagrangian in Eq. (13) 
when (j> = 1, it will be sufficient to choose < 0. With 
this choice, the kinetic term for the scalar field has the 
wrong sign but this docs not correspond to a ghost, since 
Eq. (13) is only an effective scalar-tensor action. The 
full action (11) is not plagued by any such instability at 
least in the scalar sector and at least on Minkowski back- 
ground, because there are no physical perturbative poles 
for nontrivial nonlinear interactions (e.g., [23, 25]), while 
there is a real pole for a quadratic potential.^ Notice also 



^ For a general background, the equation of motion is 

gr* [□+/(«, x)]|^ _ y/ _ leading to a very complicated prop- 
agator in momentum space which cannot be written explicitly 
except for special backgrounds. However, for a nontrivial in- 



teraction (e.g., 



n > 3) the Green equation e G{k) = 1 
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that, looking at the weak-nonlocahty Umit (13), one has We will make use of the operator identity [68] 



a = 



/MpiMV 



(14) 



where Mpi is the reduced Planck mass. Most of the stable 
solutions we shall find have \a\ < 1,^ which might indi- 
cate a trans-Planckian problem similar to that in stan- 
dard inflation, where the classical field acquires large val- 
ues above the Planck scale. However, nonlocal effects, if 
they exist, should not be negligible in the very early uni- 
verse and the approximation (13) might not be reliable is 
such regimes. Consequently, the problem stemming from 
Eq. (14) might be an artifact of the weak-nonlocality ap- 
proximation. 

Although the motivations underlying our ansatz (11) 
differ from those of [25], we eventually got a very sim- 
ilar result. This is not completely unexpected because 
both the diffusion equation and the ultraviolet proper- 
ties invoked in the construction of [25] (absence of ghosts 
and asymptotic freedom) rely on Eq. (1). The diffusion 
equation encodes some of the gauge symmetries of SET 
at spacetime level [43] and, as such, gives rise to rather 
rigid physical properties. 



B. Equations of motion 



To find the equations of motion we need tlie variations 
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= -^9^..V^Sg''\ (15) 
SR = (i?^.+5M.n- V^V.)^ff''^ (16) 



where V^V^ = di^Vfj, — ^'j^^Va is the covariant derivative 
of a vector and the curved d'Alembertian on a scalar 

6 is 



(17) 



docs not have nontrivial zeroes for any operator A and one 
has no perturbativc excitations on any background and for any 
function /. When the cigcnstatcs of A arc a complete basis 
for the functions G(k), tlie statement can be proven explic- 
itly case by case. For spacetimes with constant Ricci curva- 
ture or which arc homogeneous, / does not depend on spatial 
coordinates. Then, for a nontrivial interaction the propagator 
is G{k^) = e~^*^-^~'' ^ which, as announced, has no nontriv- 
ial zeroes. On the other hand, for a quadratic potential the 



propagator reads G{k^) = [e 



-.(/-fc^) 



-(2m2)-i + 



'^n=-oal^*(f ~ "log — '2'imi]~^ , where we expanded 
d la Mittag— LefHer. There axe infinitely many perturbative states 
with complex-mass poles — fe^ = logm^ — / -I- 2Trnr^^i. The 
theory can be rendered stable by projecting out these modes 
according to certain mathematical prescriptions [66, 67]. 
^ The only exception is Eq. (48) for p near or 2/D. However, 
solutions with these fine-tuned p's do not inflate. 



f dse'^{6X)e 
Jo 



ir,-s)X 



(18) 



for a (differential) operator X. 

Also, for two scalars $i and $2 (— ^ indicates integra- 
tion by parts), 



'-9 



-(a^$i)(9,$2)] , 



(19) 



and 



SfjR) 



'53'l*2/'(i?) 



5R 



y^{$i$2/'(i?)i?^. 

+(i?M-n-V;,v.)[$ici>2/'(i?)]}, 

(20) 

where we have discarded boundary terms (this can be 
done compatibly with the variational principle after 

adding a boundary piece to the total action). 
The scalar equation of motion 5S/S(f> = is 



</.(l-n,a;) = F'[0(l,x)]. 



(21) 



In the absence of extra matter, the Einstein equations 
are (coordinate dependence of the fields implicit) 



= 



6S 



-9 Sg^" 



= -9^.{^cl>{m^-r*)-V[cl>il)] 



= [ ds [ d^x^/^(t){l - s) 
—9 Jo J 

- r* -h s) 



5n ^6f{R) 



1 



Jo 



(22) 



^ni^is) = -g^^[(j){l - s)0(l){l - + s) 

+d„(p{l - 3)0" (pil -r^+ s)] 
-d^4>{l - s)d^4>{l -r^ + s) 

- .s)(6(l - r, + s)f{R)R^, 

+(v.n-v^v^)[(/)(i-s) 

x<^(l-r, +s)/'(i?)]. 



(23) 
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Taking the trace, 

= -D^^^<p{l)<l>{l-n)-V[cf>{l)]^+ £ dsE{s), 

(24) 

S(s) = y (/.(I - s)n0(l - + s) 
D 



^ V 2" ~ 7 a<,<^(l - s)d''(f){l -n + s) 

+cf>{l - s)<l>{l - r» + s)f'{R)R 

+iD - l)n[0(l - - n + . (25) 



III. COSMOLOGICAL SOLUTIONS 

A natural background whereon to study the model 
is Friedmann-Robertson- Walker (FRW), defined by the 
line element 

ds^ = g/^.^dx^dx" = -dt^ + a{tfgijdx'dx^ , (26) 
where t is synchronous time, a{t) is the scale factor, and 

-Qijdx'dx^ = :— ^ + g''di}l_2 (27) 

is the line element of the maximally symmetric {D — 1)- 
dimensional space S of constant sectional curvature K 
(equal to —1 for an open universe, for a flat universe, 
and +1 for a closed universe with radius a). 



A. Friedmann equations 

On an FRW background, 



Roo = -iD-l){H^ + H), 



^ + {D-1)H^ + H 



where 



R={D-1){^ + DH^ + 2H], 



H = 



(28) 
9ii, (29) 

(30) 
(31) 



is the Hubble parameter. The 00 component of Einstein's 
equations is 

= i0(l)</,(l - n) - y[(/.(l)] + J'dsEoois), 

(32) 

^oo{s) = -i(/>(l-s)n<^(l-r*+s) 

_i^(l_s)^(l-r* + s) 

+(/.(! -s)0(l-r,+ 5)/' (i?)i?oo 

+{D - l)Hdt[H^ - s)<^(l -n+ s)f'iR)]. 

(33) 



We specialize to the case 

fiR) = aR, 



(34) 



where a is a dimensionless constant. In the Friedmann 
equations (32) and (24) one has 

^oois) = ~(f>{l - s)ncl>{l - n + s) 

-l<f,{l - s)4>{l - + s) 

-a{D - \){H^ + ij)0(l - s)0(l - + s) 
+a{D - l)Hdt[(l){l - s)0(l - n + s)] , (35) 



and 



^{s) = y 0(1 - s)n(/.(l - r, + s) 



D 



-lU(l-s)<^(l-n + s) 



+Q!0(1 - s)0(l - + s)R 
+a{D-l)n[(l){l-s)(j){l-r,+s)\. (36) 

B. "Stationary" solutions with quadratic potential 

Although the diffusion method allows one to find 
asymptotic solutions for any given background, it is not 
necessary to resort to all its machinery in order to get 
some interesting results. In particular, one can look for 
simple solutions which are trivial along the diffusion flow, 
i.e.. 



(/)(r,t)=e^V(t), 



(37) 



where /? is a constant and is a stationary solution of 
the diffusion equation.^ Since /3 only changes the nor- 
malization of the field and its potential, we can set /? = 
without loss of generality: 



(□ + Qi?)(^ = 0. 



(38) 



Later on we will actually see that most of the solutions 
do require j3 — 0. 

Even if the r dependence is trivial, Eq. (38) contains 
a class of full-fledged dynamical solutions which are non- 
trivial in time and, from Eq. (21), have 



V 



(39) 



The sign of the effective potential depends on the sign of 



W{^) = (1 - e-"«) \ . 



(40) 



The "Wick-rotatcd" diffusion equation (r — >■ ir) is a Schrodinger 
equation, so /3 is the analogue of the energy eigenvalue of sta- 
tionary solutions. 
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Unless indicated otherwise, we will ehoose the sign of r, 
so that W is bounded from below. Since the interaction is 
quadratic, we regard this case as perturbative in a quan- 
tum field theory sense, although it is still fully nonlocal. 
Notice also that, for quadratic potentials, the equation of 
motion is linear in (j) like the diffusion equation. Only in 
this case does the diffusion equation contain information 
on the structure of the potential, which explains why one 
can obtain exact solutions contrary to the typical output 
of the general method. 

The Einstein equations on stationary solutions are very 
simple. For any background, S^j, = 0: 

= {^ + 2a)g^AipDip + daipd''ip)-{l + 2a)d^ipd^ip 
+aip^R^^ - 2aipV^V^ip , (41) 

where we used Eqs. (22), (39) and □((^^) = 2(pnLp + 
2dcr'pd'^ip. After using Eq. (38), the Friedmann (00 and 
trace) equations for a homogeneous field become 







" K 




— + 









1 



2a{D-l) 



= [4a(£) -1)+D-2]{lp^ + aRip^ 



(p^ + 2H(p(f , 

(42) 
(43) 



The second equation is automatically satisfied for 

lD-2 



a = —a* = — 



AD-1 



(44) 



For the time being we assume a ^ a^, and present solu- 
tions to the Einstein's equations with nontrivial ip{t) and 
a{t) profiles. 

We start with de Sitter metric, 



a{t) =e"\ H = const. 
For a flat background K = 0, one finds 

D-1 



\{D-l)Ht 



a 



AD 



(45) 



(46) 



The scalar field rolls down its potential towards its min- 
imum, reaching it asymptotically. This feature is similar 
to what was found in [18] (where, however, the Einstein 
equations were not solved) and is typical of nonlocal mod- 
els where the cosmological friction is enhanced by non- 
local operators. This is determined by the form of the 
nonlocal equations of motion, and it happens even in this 
case where the nonlocal operators are trivialized on sta- 
tionary solutions. 

Another class of solutions is power law: 



a{t) 



(47a) 
(47b) 



Inflation (a > 0) happens when p > 1. When K = 0, the 
solution is 



a = — 



{Dp-p-lf 
4(£>-l)(£)p-2)p' 



p-1). 



(48) 



In four dimensions and for p > 1/3, g < and the field 
rolls down towards the minimum. 

We now check the classical stability of these solutions 
in synchronous time formalism.^ A background solution 
{H{t),Lp{t)) is perturbed homogeneously, 

H{t) H{t) + 6H{t) , ¥'(*) ^ 'pit) + 6(p{t) , (49) 

and the equations of motion (38), (42), and (43) are lin- 
earized (the scalar equation of motion is an identity). 
The system can be written in a matrix form: 



SX = M5X, 



where 



5X = 



6H\ 
6ip ' 



(50) 



(51) 



and the entries rrig = (M)g of the 2x2 matrix M are cal- 
culated on the background solution. The characteristic 
equation 



det(M - Al) = 



(52) 



determines the eigenvalues A and a solution is stable pro- 
vided Re(A) < 0. In general, the eigenvalues are time de- 
pendent, in which case they are interpreted as evaluated 
at a given time t [71]. 

For all the solutions above mi2 = 0, so that Ai = nin 
and A2 = m22- For the de Sitter solution (46), 



Al = 0, 



Ao 



1 



{D-1)H<0, (53) 



thus implying stability. For the power-law solution (48) 

{t > 0), 



Ai = - 



t 



., = _(^^zik±i. (54) 



2t 



the latter being negative when 



P> 



D-1' 

which is true for an expanding universe {p> 0). 



(55) 



^ For an introduction to phase space analysis, see [69]. Here it is 
worth mentioning a caveat about the choice of clocks. In general, 
time t is an unphysical parameter and one has to choose an 
internal physical clock. This can be one of the matter fields in 
the total action, e.g., a scalar field or a barotropic fluid. Failure 
to do so can sometimes lead to inconsistencies in the stability 
analysis. The reason [70] is that perturbations in synchronous 
time discriminate between trajectories differing only by a shift 
in time, which are actually physically equivalent and should be 
identified. On the other hand, in the presence of an internal 
clock time shifts can be physically distinguished. In the analysis 
below, this turns out not to be an issue. One can convince oneself 
by noticing that an extra dust component in the system would 
leave the perturbed equations unchanged. 



7 



Before moving on we mention an open interpretational 
issue. To check the stability in a standard way, in princi- 
ple one should perturb the full equations of motion and 
not their push-forward on stationary solutions, unless the 
ansatz itself [Eq. (38)] is stable. Upon perturbation, the 
full diffusion equation stays linear in the scalar field but 
contains a term proportional to drScp, which we have ig- 
nored, and a source term coming from the variation of 
the metric. The problem is that, while in local gravity 
the perturbed equations are linearized and so are rela- 
tively easy to solve, here they are indeed linear, but still 
nonlocal (nontrivial r dependence). 

More work might be required in the future to clar- 
ify the stability analysis of dynamical solutions. For the 
moment, however, notice the following. Assume that a 
given solution <p[r, t) of the diffusion equation, stationary 
or not, is a solution (j){t) = (/>(ri , t) of the equations of mo- 
tion for some r = ri. In the (r, 0) plane, (/)(t) is a straight 
line at r = ri [eventually a half-line or a segment, if the 
field (p(t) is bounded]. A stationary solution has no r 
dependence, which we can fix by setting ri = 0. Suppose 
now that the linear perturbation S(f){t) = S(f){r2,t) solves 
the linearized equations for some r2. This is represented 
by another line (or segment) parallel to the first and at 
distance Ar = |r2 — ri]. For a generic nonstationary lin- 
ear perturbation, Ar will be determined by the specific 
background and perturbed solutions, and is not neces- 
sarily small. On the other hand, for a stationary linear 
perturbation Ar = 0. 

According to one interpretation, linear perturbations 
are nonstationary for a background solution with general 
self- interaction V{(j)), but in the case of stationary solu- 
tions we expect that also the perturbation is stationary 
(drS(f> = 0). The main reason is that the diffusion equa- 
tion is an external input which knows about the dynam- 
ics only via the background metric, but that in general 
is ignorant about the shape of the scalar potential. How- 
ever, for stationary solutions the background potential V 
is quadratic, which is the only case where the diffusion 
equation is implicitly sensitive to the scalar field dynam- 
ics. In fact, the diffusion equation can encode informa- 
tion about the potential by adding a mass term, but this 
is a trivial normalization along the diffusion flow. There- 
fore, the background and perturbed diffusion equations 
have basically the same implicit knowledge of the scalar 
dynamics. The theory of diffusion suggests that the only 
exact solution S(f) will be stationary as well or quasis- 
tationary (Ar ^ 0), due to the delicate dependence of 
the solution of the diffusion equation from the boundary 
conditions. 

A second interpretation defines the stability as per- 
taining only to the spacetime dependence and not, as in 
nonstationary perturbations, also to the r dependence. 
This is consistent with the idea that the r dependence 
of the scalar is really a mathematical trick, and it is the 
same as imposing that the ansatz (38) be stable by defini- 
tion. In the case of nonstationary backgroimd solutions, 
this option imposes that the most general perturbation 



of the most general solution (f){ri , t) be always of the form 

The second interpretation prescribes that all solutions 
be perturbed by stationary fluctuations, while the first 
reaches this conclusion only for stationary solutions and 
only at the linear level. In the above stability analysis 
we have considered stationary perturbations of the form 
(49) and tacitly embraced the second view. 



C. General solutions with quadratic potential in 
conformal gravity 

The flat de Sitter solution can be generalized straight- 
forwardly to any f{R) theory. In fact, it is suflicient to 
make the replacement 



a aca = /'(-Rds) 



(56) 



in Eq. (46), where R^s = D{D — 1)H^. Depending on the 
form of /, the value of the Hubble constant is determined 

by ckoff. 

Another generalization is based on the fact that the 
most general solution with a = — a* can be found for 
any FRW background via the following shortcut. 

The critical value Eq. (44) is well-known in conformal 
gravity models [72]. In four dimensions, a = —1/6. Con- 
sider a metric g^i, and the conformal transformation 



for some r-independent fl 
bols transform as 



(57) 

n^'^cl), (58) 
--n{x). The Christoffel sym- 



■ IJ,!/ 



(59) 



where 



= V^lnO, 



= V''Ou = ninfi. 



(60) 



Noting that 

Q.-'^R = R + 2{D-l)0-{D-l){D-2)0^0^' 

Q.-^U = □ - (D - 2)C'^V^ , 



(61) 
(62) 



one can see that the combination (□ + aR)4> has a well- 
defined conformal weight only if a = — a*: 

(□ - a*R)(t) = fi^(Q - a».R)0 . 

i?o = const, the diffusion 



On a background where R 
equation is equivalent to 



(63) 



We specialize now to FRW backgrounds, where the met- 
ric can be written in conformal time 



dt 



(64) 
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so that n = fl-i, □ = -52, and Rq = 2{D - 1)k. For 
stationary solutions, the diffusion equation reduces to 
Eq. (63): 



if" + sgn(K)7V = 



where 
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D 
2" 



f 



(65) 



(66) 



and primes denote derivatives with respect to t. Just 
solving this equation, one obtains a very wide class of 
solutions from Eq. (58), 



(67) 



In fact, the whole dynamics reduces to one equation, en- 
coding both Eq. (65) and a constraint on (p' and the 
curvature. To show this, instead of the trace equation 
(43) one takes the ii-component equation Yin = 0, 



2a 



iRip'^) +aip'^R + 2aHipip. (68) 



Replacing Eq. (67) in Eq. (42) or (68) with a = —a*, one 
obtains just the same equation 

(^')'+sgn(K)72^2 ^0, (69) 

which can be solved only ifK = OorK = — 1. There 
are no closed-universe real solutions. Differentiation of 
Eq. (69) in conformal time yields Eq. (65), so for the 
critical value (44) the stationary problem is drastically 
simplified. This is possible only for stationary solutions, 
where geometry is factorized out of the equations of mo- 
tion. Let us consider a few concrete examples. 



1. Flat universe (k = 0^ 

The g frame is Minkowski and = 1 is a solution of 
the free equation if" ~ 0.^ Then, the general solution is 

^(r) = [a(r)]i-^. (70) 
For the de Sitter background (45) in synchronous time, 

<p{t) ^ e-i^-')"\ (71) 

and the evolution of the scalar field is qualitatively the 
same as for Eq. (46). 

For the power-law profile (47) , the solution has 



1 p. 



(72) 



^ Also ip = T solves the diffusion equation, but it does not solve 
Eq. (69). 




FIG. 1: Bouncing K = solution (74) with D = 4 and Ho = 2. 
Thin line: scale factor a{t). Thick line: scalar profile ip{t) = 
l/a{t). 



Since the algorithm given by Eqs. (65) and (67) is valid 
for any choice of the scale factor, there is an infinite set of 
solutions without big bang singularity. A D = 4 bounc- 
ing solution with zero intrinsic curvature is (r* > 0) 



a{t) 



cosh t , 



i? = 6[1 + (tanht)^] 



H = tanh t . 



(73) 



In the contracting phase the scalar field rolls from its 
minimum up to a maximum value. At the inversion point 
the universe bounces and starts expanding, while ip roll 
back towards the global minimum. This solution is the 
scalar-tensor analogue, in the sense of Eq. (11), of the 
bouncing solution of [25]. The fact that the scale factor 
a is nonsingular, however, is to be ascribed more to the 
choice of a conformal operator rather than to the good 
ultraviolet properties of these nonlocal theories. 

Another bouncing profile is the superaccelerating cos- 
mology 



a{t) = exp -Hot 



1 



H — Hot , 



(74) 



which is plotted in Fig. 1 for 13 = 4. The scalar field 
evolves from the global minimum up to the inversion 
point ip — 1 at the bounce. Then, cosmological friction 
drags it back to the minimum. 



2. Open universe (k — —1) 

In the open case, Eq. (65) is solved by exponentials and 
there are two general dynamical solutions which read 



V'±(t) = [a(T)] ^ e 



7 pT7T 



(75) 



-1.0 -0.5 0.0 0.5 1.0 




-2 2 4 

t 



FIG. 2: de Sitter K = -1 solutions <^±(f) [Eq. (76)] with 
D = 4. 



Therefore, and somehow surprisingly, there are two "de 
Sitter" [H = const) solutions also for an open universe: 

<P±(t) = e-^'^*exp(±-le-^*) , (76) 

where we used r = —e^^*/H. The profile has the 
same features as the flat solutions and rolls towards the 
global minimum. The solution Lp- climbs from the min- 
imum up to some maximum value, then rolls back down 
(Fig. 2). 

For the power-law profiles (47), when p = 1 (linear 
scale factor) the solutions have 

In four dimensions, only the positive root is nontrivial, 
q = —2. For arbitrary p, the solutions are 

<^±(i) = t-^'Pe^^*'"'. (78) 

The typical plot of, say, ip- is shown in Fig. 3 and is sim- 
ilar to the previous one. For comparison, the exact de 
Sitter and power-law solutions in standard general rela- 
tivity are reported in the Appendix. 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

t 

FIG. 3: Power-law K = -1 solution if -it) [Eq. (78)] with 
D = 4 and p = 3. 
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2 - 



1 - 



0^ 




-2-10 1 2 

FIG. 4: Asymmetric bouncing solution given by Eqs. (75) and 
(79), with D = 4 and Ho = 2. Thin line: scale factor a{t). 
Thick line: scalar profile ip-(t). 



For the bouncing solution Eq. (74), the conformal time 

is 




thus leading to the asymmetric lump [t) shown in 
Fig. 4. The evolution is similar to the flat case, except 
that now open geometry helps the scalar field to climb up 
to a higher inversion point before being dragged back by 
cosmic superfriction. The lump ip+it) has its maximum 
at f < 0. This corresponds to a solution with a different 
initial condition for the scalar field: ip+ rapidly reaches 
the (higher) inversion point before the bounce. 

Still in an open universe, we sketch some other solu- 
tions in four dimensions which well illustrate the exotic 
properties which can emerge in a nonlocal conformal set- 
ting. The first case is periodic in time, so we restrict it 
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(a, 1^) are fixed points and there are no attractors. This 
is a rather bizarre situation in cosmology. Typically, a 
solution of the background dynamics is probed by per- 
turbing it slightly in all sectors, inclusive gravity, as done 
previously. Thus one checks whether one will hit the 
same solution after evolving the system from slightly dif- 
ferent initial conditions. In the conformal system with 
a = —a*, on the other hand, this operation is ill defined 
because a{t) is a fixed input. 

Can we conclude that the scale factor is nondynami- 
cal and, thus, the conformal solutions of this section are 
not physically sensible? We argue in favor of a negative 
answer. Nonlocal models are notoriously rigid and the 
choice of initial conditions is much more restricted with 
respect to local theories with a standard Cauchy problem 
[4, 10, 73]. Here we see this property in action, with a fur- 
ther restriction of the background choice. Moreover, the 
scale factor factorizes out only because we have assumed 
trivial diffusion, Eq. (65). A realistic stability analysis 
should go beyond the homogeneous level, where the so- 
lutions are no longer stationary in the diffusion sense; 
this would correspond to looking at semiclassical (inho- 
mogeneous) perturbations. 



D. Asymptotic solutions and spontaneous 
symmetry breaking 



FIG. 5: Upper panel: scale factor lOa(t) (thin line) and scalar 
profile ip{t) (thick line) for the recoUapsing solution Eq. (80). 
Bottom panel: scale factor a{t) (thin line) and scalar profile 
ip{t) (thick line) for the expanding solution Eq. (81). 



to half a period with positive scale factor {r^, < 0): 



a{t) = sint . 

K = -1, 



t e [0, tt] 

i? = -12, 



(80) 



1 



1 — cos t 



which is shown in the upper panel of Fig. 5. Despite being 
open, the universe recollapses onto itself while the scalar 
field rolls towards the global minimum ip — 0. Eventually 
one hits a future big crunch singularity at finite times. 
Another solution consists in a "Wick rotation" of the 
former (r^ > 0): 



a{t) — sinht , 
K = -l, R = 12, 
1 



(81) 



cosh ^ — 1 ' 



depicted in the bottom panel of Fig. 5. 

The usual stability analysis does not apply to the above 
solutions because a{t) is factorized out of the equations 
of motion. Thus, all points in the phase space plane 



We have found exact, cosmologically nontrivial solu- 
tions which are constant along the diffusion direction but 
one might ask whether there exist also nonstationary dif- 
fusing profiles. This is a rather difficult question to an- 
swer because there is no systematic method to solve the 
equations of motion starting from a solution of the dif- 
fusion equation. In this paper we will limit ourselves to 
some general remarks for f{R) = aR. 

The diffusion equation transfers the nonlocal degrees 
of freedom into the extra direction [40, 43] and one can 
study the system locally in spacetime coordinates with- 
out risk of falling into the contradictory situations typical 
of nonlocal dynamics. In particular, it is possible to ex- 
pand field profiles and equations in series and look for 
asymptotic solutions. 

The static potential felt by the scalar particle is, ac- 
cording to Eq. (11), 



(82) 



The curvature term matches a possible quadratic term 
m?4>'^ 12 in V and the net result is a dynamical rescal- 
ing of the squared mass, which may be even negative 
(tachyon) . A particularly important situation is when V 
has a minimum at the origin with a power higher than 2 
(e.g., V (fy^ but it could be of any other form). In this 
case, the dynamical term of the potential always converts 
the minimum (j) = Q oi V into a local maximum of W ^ 
and the minima are located at 0o 7^ 0. Classically, stable 
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solutions should tend to these minima, 



(83) 



Quantum mechanically, the scalar field is expected to un- 
dergo spontaneous symmetry breaking and take a non- 
vanishing expectation value on the vacuum, {(j))o = (f)Q. 
In this case the system relaxes to pure gravity, Eq. (12), 
the expectation value of the scalar field playing the role 
of Newton's constant. This is precisely the same kind of 
mechanism that happens in string theory with the dila- 
ton field. 

If the scalar field potential is to be determined a pos- 
teriori, in the presence of gravity a convenient approach 
is to factorize V out and consider the equation of motion 



0: 



/ * rfs(DEoo 
Jo 



(84) 



which can be solved for given FRW profiles. In turn, this 
determines the scalar field profile in the diffusion equa- 
tion. Looking for asymptotically (nonvanishing) con- 
stant solutions for the matter field in Eq. (84), one gets 
DRoQ + R = Q, that is, 



2k 
D-2 



(85) 



This constrains the asymptotic geometry of spacetime. 
When K = 0, the only solution is de Sitter. For K = ±1 
and D = 4 the scale factor is, respectively, a = (cosh bt)/b 
and a = (sinh &<:)/&, where b is an arbitrary constant.^ 
The curvature is always a positive constant, Rq = 126^. 
In the limit 6 — >■ 0, the K = — 1 solution reduces to the 
linear case a — t. This is the only nontrivial solution 
of Eq. (85) with asymptotically vanishing curvature i?o- 
Whatever the full solutions are, they can only be con- 
structed by considering the above scale factors as either 
exact or asymptotic profiles. 

We conclude this part by determining the relations be- 
tween the constants defining the spacetime geometry (A 
and G) and those coming from matter [Vq = V{(j)Q) and 
(/)o]. This can be achieved in two independent ways, by 
looking either at the action Eq. (11) or at the equations 
of motion, Eqs. (22) and (23). Taking the weak non- 
local limit of the action (r, ^ 1), we get the standard 
Einstein-Hilbert action with cosmological and Newton 
constants given by 



A, 



2r*a 



21^ 



87rGw 



(86) 



By consistency, we expect Eqs. (22) and (23) to repro- 
duce the Einstein equations in vacuum with the values 



(86). To this aim, it is important to note that the asymp- 
totic value of the scalar field depends on r* . The potential 
W is explicitly r* dependent and, as a consequence, so 
are its local minima. 

Once again, we can appeal to the diffusion equation 
governing the flow in the r* variable. In fact, 

0(r-r*,a;) = (j){r,x) — rt,dr(j){r,x) + 0{r1) 

= [l+r^U + r,aR + 0{rl)](t){r,x){B,7) 

One could have just ignored the difFiision picture and 
considered a small r* expansion in 0(1 — r*,a;) = 

(l){l-r^,x) = [l + r^n + r^aR + 0{rl,n'^)\<i){l,x) . (88) 

However, in the latter expression it is not clear whether 
we are entitled to safely take an asymptotic limit 0asymp 
for (j), as we have to know the contribution of all the 
derivatives of the field. On the other hand, the diffu- 
sion method regards 0(r, x) as a local field with two ar- 
guments, and the small r* expansion is a genuine local 
expansion. In other words, the diffusion picture justi- 
fies manipulations of asymptotic solutions in the nonlocal 
model. 

Taking the asymptotic limit a; ^ oo of Eq. (87), we 
get 



lim ( 



n,x) = [I + r^aRn + 0{rl)](t)f, 



(89) 



Performing the asymptotic limit in the equations of mo- 
tion (22) and (23), expanding in r* up to the first or- 
der, and taking Eq. (89) into account, one easily gets the 
standard Einstein equations in vacuum, with Newton and 
cosmological constants as in Eq. (86). 

Another interesting avenue to explore for nonstation- 
ary cosmological solutions is the conformal case a = — a*. 
Then, the diffusion equation becomes 



(j)" + sgn(K)7^(^ - a^~^dr 



0. 



(90) 



Interestingly, in six dimensions this is the usual homo- 
geneous diffusion equation in conformal time, of which 
we know the solutions. For example, on a de Sitter back- 
ground the Gaussian profile (j){r, t) becomes an asymmet- 
ric kink in synchronous time, while for the power-law ex- 
pansion it is a very flat lump. One should check, however, 
if the profiles (j> solve, even asymptotically, the Einstein 
equations for a given potential. This is not guaranteed. 
In fact, unlike the case of stationary solutions the system 
is not trivialized to Minkowski. The reason is simple. Let 



hn = (□-a*i?0"<^- 



(91) 



The field 0„ is conformal with weight a~^" when n = 
(D/2 — l)/2. Hence, none of the is conformal (except 
possibly one if D = 6 + 2k, fc e N) and the object 



Clearly, the symmetry breaking argument does not apply to the 
corresponding exact stationary solutions with quadratic poten- 
tial, Eqs. (46), (73) and (81). 



(92) 



n=0 



contains an infinite number of a factors, which survive in 
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E. Inclusion of cosmological fluids 

So far we have been focusing on solutions where the 
only matter component was the nonminimally coupled 
scalar field. The first reason to do so was simplicity; the 
structure of the above solutions is already rich enough 
in the stationary "empty" (no fluid) regime to illustrate 
some features of the class of nonlocal theories under in- 
spection. Second, as mentioned in the introduction, we 
expect nonlocality to manifest itself in the early uni- 
verse, in particular, through a bounce or during inflation 
when dust and radiation are typically (but not always) 
"switched off." Third, our purpose was to define a non- 
local theory which would be consistently nonlocal in all 
sectors. This means that it is not clearly consistent to 
introduce a traditional cosmological fiuid, because its ac- 
tion is intrinsically local. Instead, nonlocality can be 
easily and rigorously implemented in a field theory and, 
if one were interested in building a realistic multicompo- 
nent postinflationary cosmology, one should define a set 
of fields with asymptotically perfect-fluid-like equations 
of state. 

Of course one can just ignore this issue, add ordinary 
fluids, and study their cosmology numerically, eventu- 
ally placing observational constraints on the parameters 
of the theory. This exercise lies beyond the above men- 
tioned foundational framework, but in this section we can 
make some general comments on analytic single-fluid so- 
lutions. 

Let us add to the action (11) a generic matter part 
S'm giving rise to the energy-momentum tensor T^i, = 
-{2/^f^)5Sra/Sg^"'. The latter constitutes the left- 
hand side of the Einstein equations (22). For an ordinary 
perfect fluid with energy density p and pressure p. the 
energy-momentum tensor is T^^ = {p +p)UnUi, + pgij,^, 
where is the comoving four- velocity field (unit timelike 
vector, Uf^u^ = —1) tangent to a fluid element's world- 
line. The continuity equation reads p+{D — l)H{p+p) = 
0, while Eqs. (32) and (24) become 

p = ^<^(1)</.(1 - n) - 1/[0(1)] + £ ds Eoo(s) , (93) 



{D-l)p-p = -D{ -<l>im'L - r*) - V[m] 



+ 



[ ' ds E(s) . 
Jo 



(94) 



On stationary solutions, the Einstein equations are 



T = 



- -h 2aJ fif^^ {^pUip + daipd'^ip) 

-{1 + 2a)dfj,ipd^(p + aip^Rn^ 
-2aipVfyvip, 



where we set r, = 1 in the trivial integration in the right- 
hand side. In particular, the Friedmann equations are 



p = a{D- 1) 

+2a{D - VjH^if 



2 -2 

^ -2^ 



p-{D-l)p = 



2a{D -1) + 



D 



(96) 
(97) 



For a barotropic fluid, the equation of state is p = wp, 
where w is constant. Then, from the continuity equation 
it follows that p = poa~''^~^^^^~^^\ For pressureless dust 
w = 0, while for radiation w = 1/{D — 1). We would 
like to see how some of the previously found solutions are 
modified in this setting. We take po > and — 1 < w < 1. 

As before, let us first assume a ^ —a*. For the de 
Sitter metric (45) and a flat background, the solution 
(46) is replaced by 



^it) = ^^e-iiD-i)i^+i)Ht ^ 
(L»- l)(l-w;2) 



a 



AD 



(98) 



where 



± 



Dpo 



{D - 1)H y -w{w + 1)[1 + {D- l)w\ ' 



D 



(99) 



Neither radiation nor dust are solutions. The scalar fleld 
rolls down the minimum of the potential but at a slower 
rate with respect to the no-fluid case, by a factor of 1 + w. 

The K = power-law solution (47) and (48) now be- 
comes 



a 



{l-w)[{D-l){w + l)p-2] 
4:{Dp - 2) 

-kD-l){w + l)p, 



(100) 



(95) with 
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±2, 



' (£'p-2)po 

[1 - (£) - l)wp][{D - l)wp + p- 2][{D -l){w + l)p - 2] ■ 



(101) 



r 



The solution is real only for a certain range of the pa- 
rameters. For instance, there are no real solutions for 
radiation, w = !/(£>— 1), while iov w = Q it must be 
either p>2or2/£><p<2/(£)-l). Notice that now, 
in four dimensions, g < if p > 2/[3(w + 1)]. In order 
for the scalar field to roll towards its minimum, and as- 
suming < w < 1, p must be larger than the value 1/3 
obtained in the case without fluid. 

The conformal case a = —a* is straightforward. The 
trace equation vanishes identically only in the presence of 
a radiation component, w = 1/{D — 1). Then, both the 
00 and ii components of the Einstein equations reduce to 



{^'Y + sgn(K)7"^" + 2po = 



(102) 



which generalizes Eq. (69). The flat background no 
longer admits solutions and only open-universe profiles 
survive. Contrary to Eq. (75), the general solution is a 
nontrivial superposition of the solutions of the diffusion 
equation: 



<p(t) = ¥'o[a(T)] 



-1 



27Vg 

(103) 

As before, the "de Sitter" {H = const) and power-law 
solutions are qualitatively similar. The scalar always rolls 
towards the global minimum, monotonically if po/Vo ^ 
and with a smooth, temporary turnaround if po/^pQ '> 1. 



IV. STRINGLIKE ACTION AND COSMOLOGY 

The model studied so far is based on a gravity exten- 
sion of a p-adic-like scalar field action, which has a trivial 
local limit r* — >■ (no dynamical degrees of freedom) . We 
can now make a slight but crucial modification of the ki- 
netic operator e'^ into a transcendental expression of the 



form e'^JC 



S = 
where 



.= p + /(i?)]0. 



(104) 



(105) 



The linear case f{R) = aR is important for several rea- 
sons already illustrated in [25]. First, when r, < one 
formally recovers the effective spacetime action of the 
tachyon and the graviton in closed SFT. When gravity is 
switched off, Eq. (104) reads 



\r,\a 



□<^-y(0) 



(106) 



where </> is (up to a (jxj) negative mass term) the tachyon 
field "dressed" with an exponential operator (see, e.g., 
[42] for details) and F is a polynomial potential. On the 
other hand, when ^ = 1 we get 



-lr.|(n+aH)^_ 2a' 



(107) 



This is an effective nonlocal gravitational action in terms 
of Riemann invariants. In closed SFT we do not have 
any such thing because the effective level-truncated ac- 
tion for Qfj^i, is obtained in a nondiffeomorphism-invariant 
form [47, 58].^ However, by expanding Eq. (107) at lin- 
ear order around Minkowski spacetime, one indeed gets 
the correct propagator G{k^) in momentum space for the 
graviton, as shown in [25]: 



(108) 



This propagator is ghost free and realizes asymptotic 

safety at large momenta, two properties expected in 
a ultraviolet-finite nonperturbative theory of quantum 
gravity [4, 23, 25].^ Therefore, Eq. (104) is interesting 
both as a toy model for the spacetime effective dynamics 
of gravity in closed SFT and as a nonperturbative ansatz 
for the gravitational action at high energies/large curva- 
ture, which can play a major role during the very early 
universe. 

From Eq. (104), the scalar equation of motion 6S/6(p = 
is 



4>{i-n,x) = V'[cl,{i,x)], 
while the Einstein equations read 



(109) 







= -5^. 1^.^(1)^(1 -n)- mi)]} 

+S^,(0)+ / *dsS^,(s), 



(110) 



where S^,^ is J^^i, with 0(1 — r* -|- s) replaced by — 
r* -|- s) (or the symmetrized expression). With obvious 



® Nonetheless, diffeomorphisms are still part of the symmetry 
group of the theory [47]. 

® This type of propagator appears also in superrcnormalizable 
quantum field theories [74] and in the context of noncommutative 
theories where a minimal length is effectively induced [75, 76]. 
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notation, the trace equation is 



= <! -0(1)^(1 - r*) - 1/[0(1)] 



m 



ds S(s) . 



(Ill) 



The equations of motion could have been obtained also 
by replacing V r*F in the equations of the p-adic- 
like case, Eqs. (21) and (22), and differentiating with re- 
spect to r*. One notices, in fact, that dr,(f>{l — r^, + s) = 
4>{1 — r^, + s) and that the contribution of the upper ex- 
trcmum of the integral, S^j,y(r*), is equal to 2^,^(0) upon 
symmctrization (i.e., splitting JS in two and replacing 
s — )■ r, — s in one of the pieces) . 

To get stationary solutions with f{R) = aR, it is nec- 
essary to reinstate the normalization constant 13 in the 
diffusion equation, 



^ + 13^ = 0, 
which is fixed by Eq. (109): 



(112) 



In the presence of a barotropic fluid, there is a richer 
structure of exact solutions. The equations of motion 
read 



p = a{D - 1) 
1 



2a{D - 1) 



—iP" + 2Ha{D-l)ipip, 



(120) 



2[1 - (£> - l)w]p = [Aa{D - 1) + D - 2]{ip^ + aRip^) 
+ [4a(£'-l) + £>]/3^^ (121) 

The scalar field profile of the de Sitter solution now is 
<^(i) = <^oe-^(^-i)('"+i)^*, (122) 



with 



a 



(^-1)^(^ + 1) + 7?^ 
' A{D- l)[D+{D- l)w] ' 

■-(D - l)[ADa +{D- 1)(1 - w;2)]ir2 _ 



(123) 



2 ^ 



(113) 



The equations of motion are still S^i/ = but with ^ 
(3^1. The Priedmann equations read 



= 



+ 



y ~ V 2a{D-l) 

[Aa{D - 1) + i:* - 2] (<^2 _^ aRip^) 

+ [4a(D - 1) + D]l3if^ . (115) 



The requirement /3 ^ is very stringent, a ^ —a*, and 
the only exact solution in common with the "j)-adic" case 
is the flat (k = 0) de Sitter proflle (45) and (46) with 

p ^ a(4Z^n +D-1) [MP - 1) + ^ - 2] 



(4a +1)2 



^p{t) = exp 



2aH 
4a + 1 



-t 



(116) 
(117) 



The solution is stable only for certain values of a. In 
fact, the eigenvalues of the characteristic equation for 
the linearized system (50) are 



Ai = 



1 



D \ H . 



,4a + 1 

The solution is stable for 
D-l 



A, 



2aH 
4a + 1 



AD 



< a < 0. 



(118) 



(119) 



There are no constraints on (/' since a diverges at the 
exotic value w = —D/{D — 1) < —1. 

The extra fluid terms compensate those proportional 
to /3 and allow us to consider also the conformal case 
a = — a*. The equations of motion are solved for 



[1 - (£> - l)w]po 



P = 



(124) 



where the scale factor is 

e±^t + A{D - 2)K(^^e=F^* 



a±{t) = 
A = 



4^2[-l - {D - l)w]po ' 

Ik, 2 



-1 -(£>-!)«; ifo' 



(125) 



Notice that w < — 1/(_D — 1) and neither dust nor ra- 
diation are real solutions. When K = 0, the solution 
a± reduces to de Sitter. Setting (p^ = l/[4(£> — 2)], the 
other two cases are easily written. The open-universe so- 
lution is a a siuhAt and expands monotonically, while 
for K = 1 there is a bounce, a oc cosh At. For other values 
of (po the bounce is asymmetric. 



V. DISCUSSION 

In this paper we have constructed and solved, on cos- 
mological backgrounds, an effective nonlocal model of 

gravity nonminimally coupled with a scalar field. The 
actions (11) and (104) are nonperturbative both in the 
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order of curvature invariants/'^ 

e^~l + i?+^i?2 + ... , (126) 

and in the number of derivatives acting on the metric, 

e^Rr-^ R + nR+^a^R+ ... . (127) 

Nonlocal cosmology is radically different from higher- 
order cosmological models of f{R), Gauss-Bonnet, or 
/(Gauss-Bonnet) gravity. This is because truncation of 
a nonlocal model in spacetime derivatives produces an or- 
der n Ostrogradski-like problem with altogether different 
physical properties. The well-known fact that theories 
with an infinite number of derivatives are not the large n 
limit of finite-order actions has been also invoked to ques- 
tion the relevance of higher-order cosmologies in the early 
universe [79]. Rather than a finite-order truncation of the 
gravitational action, near the big bang curvature effects 
should be consistently taken into account only within a 
fully nonperturbative framework in the sense of the left- 
hand side of Eqs. (126) and (127). The right-hand side 
of Eq. (127) might not be even well defined on a gen- 
eral nonlocal solution of the system [10]. The diffusion 
equation method allows one to deal with the full nonlocal 
operators and bypass the problems of a series expansion. 

The diffusion structure we have explored is asymmetric 
in the gravity and scalar sector; in fact the former does 
not diffuse at all. The only solutions we have been able 
to find do not diffuse even in the matter sector (more 
precisely, they are stationary along the diffusion flow), 
but in general a nonstationary diffusion structure is nec- 
essary to solve the system with a self-interacting (higher- 
order potential) scalar field. For the purpose of finding 
analytic solutions, this should exclude the a priori as- 
sumption that, preferring a "symmetric" formulation of 
the model, also the scalar sector does not diffuse. In 
this case, geometry through the curvature term ,f(R) 
would replace diffusion along r. Therefore, the theory 
of diffusion associated with nonlocal actions would be 
simply defined differently: Diffusion always takes place 
through geometry, but in the case of trivial geometry 
(Minkowski background), this is realized by an auxiliary 
higher-dimensional structure. If this was really the case, 
however, it would be probably difficult to find analytic 
or semianalytic solutions with nonlinear self-interaction 
(nonquadratic V). 

The actions we have studied are structurally similar to 
the one advanced in [25] for the following reason. On one 
hand, the proposal of [25] aimed at an ultraviolet-finite 
action for quantum gravity which would address the big 
bang singularity problem. On the other hand, we wanted 
an action which would be nonlocal in both matter and 



gravity sectors and be endowed with a diffusion struc- 
ture allowing one to reduce the dynamics to a set of local 
equations with both a second-order differential structure 
(in spacetime) and an algebraic structure (in the diffu- 
sion direction [10, 40, 42]). These questions, however, 
are implicitly related: ghost and asymptotic freedom are 
determined by the specific choice of pseudodifferential 
operators, in this case one with a natural diffusion struc- 
ture. So diffusion and good ultraviolet properties are tied 
together, as expected in string field theory [43]. 

There is, anyway, a caveat in this comparison. We not 
only stressed the importance of solving a fully nonlocal 
action with both gravity and matter cosmological non- 
trivial profiles, but in doing so it was also shown how 
these profiles can differ, even considerably, with respect 
to local scenarios. The simplest cosmological profiles 
(de Sitter and power law) are exact solutions of the non- 
local dynamics. There are a couple of remarkable facts 
associated with that. First, we needed only to look at 
stationary solutions along the diffusion flow. Second, 
contrary to standard general relativity these profiles cor- 
respond to exact dynamics even when the intrinsic cur- 
vature K is negative definite. In particular, de Sitter is an 
exact solution for a nonconstant scalar field profile, also 
in an open universe. When the nonlocal operators are 
chosen to be conformal, for models with f{R) = —a^,R 
we have found the general solution for any flat or open 
FRW background, embodied by Eq. (67); these results 
have been extended to the inclusion of a barotropic fluid, 
Eq. (103), also for stringlike actions. Within this class 
there are solutions without big bang singularity, but there 
also exist an infinite number of solutions with big bang. 
Therefore we incline not to link nonsingiilar solutions 
with the ultraviolet structure of the nonlocal action. 

At any rate, the space of solutions is likely to be much 
larger than the portion we have explored here. All our 
exact solutions have a quadratic potential. Highly non- 
linear equations of motion are of great interest, especially 
in string theory, but the exact solutions can give some 
indication of the behaviour for general potentials near 
a local minimum. Cosmological friction modifies the dy- 
namics of nonlocal scalars with respect to Minkowski and, 
in particular, should drastically change the rolling of the 
tachyon in string field theory. 



Appendix: de Sitter and power-law solutions in 
standard local cosmology 

Consider a £>- dimensional universe filled only with a 
scalar field with potential V and A = 0. The standard 



Actions of the form (126) are sometimes dubbed "exponential 
gravity" [77, 78]. 



Deviations from local cosmology is not limited to background so- 
lutions. It would be interesting to study the inflationary spectra 
stemming from the inhomogeneous perturbation of the Einstein 
equations. 
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Friedmann and continuity equations are 

D-l\ 2 



V 



,2 ' 



+ H 



D~ 1 



D 



-J -r 



(A.l) 

(A.2) 
(A.3) 



In de Sitter, Eq. (45), for a flat universe the exact solution 
is just a cosmological constant. 



m = 00 , v{<i>) = 



(£)- l)(D-2)i?2 



K = 0. 



The Friedmann equations show that there is no solution 
if K = — 1, while there is one for a closed universe, but 
only in Z) = 4: 



-Ht 



K = 1. 



This solution is not de Sitter because spatial sections 
are not flat. (We recafl that "de Sitter" in a cosmolog- 
ical sense, H = const, corresponds to the mathematical 
de Sitter spacetime only if spatial sections are flat. In 
that case, the de Sitter hyperboloid is only half covered 
by FRW coordinates.) The continuity equation fixes the 
potential: 



vi<t>) 



H 



2 j,2 



(A.5) 



The scalar field (j)± rolls down its potential from t = —oo 
and climbs it again after passing the global minimum. 
The solution is actually unique, since cosmological equa- 
tions of motion are invariant under time reversal, and it 
does not matter the direction of the rolling in a symmet- 
ric potential. 

For a power-law expansion. 



a{t) = tP , 



(A.6) 



one can try the profile (f){t) = {(j)o/q)t'^ in the (sum of the) 
Friedmann equations, but one soon finds that it must be 
q — 0. This suggests to consider the limit 5 — )• 0, which 
is a logarithmic profile: 



This gives 



(l>{t) = (/-olnt. 



{D-l)p-l o_2^/^„ 



(A.7) 



(A.8) 



(A.4) If the universe is flat. 



= ± 



{D - 2)p 



K = 0, 



(A.9) 



while for a curved universe the case p = 1 is the only 
solution: 



= ± 



D 



2k 



1. 



(A.10) 



This solution is real if D > 2(1 — k). Therefore, it is 
always valid for a closed universe, while for an open uni- 
verse it exists only in D > 4. 

In the great majority of applications in the litera- 
ture, the curvature is ignored because its contribution 
is washed away by inflation. However, it is interesting 
to note that the only curved solution in "de Sitter" is a 
closed four-dimensional universe, while for a power-law 
expansion both signs of the curvature are allowed but for 
D > 5. In D = 4, only the closed solution is allowed. 

The nonlocal cosmologies described in the main body 
of the paper, on the other hand, only allow flat and open 
solutions, but without constraints on the dimensional- 
ity of spacetime. In this sense, nonlocal cosmologies are 
"complementary" to the usual ones. 
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